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Abstract 

In this paper we study the deformations of bihamiltonian PDEs of hydrodynamic type with 
one dependent variable. The reason we study such deformations is that the deformed systems 
maintain an infinite number of commuting integrals of motion up to a certain order in the 
deformation parameter. This fact suggests that these systems could have, at least for small 
times, multi-solitons solutions. Our numerical experiments confirm this hypothesis. 
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1 Introduction 



The main purpose of this paper is to study the effects of a "deformation" of an infinite dimensional 
completely integrable system. 

The first point is to define what a deformation is. In this paper we deal with bihamiltonian 
systems of hydrodynamic type for which it is quite natural to define deformations in terms of the 
Jacobi identity: the deformed bihamiltonian structure satisfies the Jacobi identity only up to a 
certain order in the deformation parameter. 

The interesting deformations are the deformations that cannot be obtained from the original 
bihamiltonian structure by just a change of coordinates. Therefore the first problem is to select 
the non trivial deformations. 

In order to solve this problem it is convenient formulate it in terms of Poisson cohomology. 
Recently Degiovanni, Magri and Sciacca (see [4]) proved that the first two Poisson cohomology 
groups of a Poisson manifold (M, P) are trivial when M is the loop space {S 1 — > M™} and P is 
a Poisson bracket of hydrodynamic type. Getzler independently (see [11]) proved that all groups 
H l (P, M) for % positive are trivial for such (M, P). 

This result, as we will see, simplifies remarkably our problem and allows us to solve it (in this 
paper we classify the deformations up to fourth order). 

For the second order deformations we show explicitly how to obtain an infinite "hierarchy" of 
hamiltonian equations . We will see that the corresponding flows commute up to the order of the 
deformation. 

One typical class of solutions of infinite dimensional completely integrable systems is the class 
of multi-soliton solutions. A natural question arises: do the equations of the deformed hierarchy 
have multi-solitons solutions (at least for small times)? 

The numerical experiments we have performed for an equation of the deformed hierarchy show 
the existence of solutions analogous to two-solitons solutions. 

Finally we observe that the deformations of the bihamiltonian structures of hydrodynamic type 
appear in the framework of Frobenius manifolds where, with some additional constraints, they 
play a crucial role in the problem of reconstruction of a 2D TFT from a given Frobenius manifold 
studied by Dubrovin and Zhang (see [8]). One of these constraints, called quasi-triviality, is anal- 
ysed in the last part of this paper. 

The paper is organized as follows. 

The first part (section 2) is a brief introduction to Poisson cohomology. We focus our attention, 
in particular, on the infinite-dimensional version of Poisson cohomology. To do this we use the 
formalism of formal calculus of variations (see for example [10]). The main purpose of section 2 is 
to explain how to get the formulae for the Schouten brackets that will appear in the calculations. 
In section 3 we give the classification of deformations up to fourth order in the deformation 
parameter. 

In section 4 we construct the deformed hierarchy and we find one soliton solutions of one equation 
of the hierarchy. 

Section 5 gives the results of the numerical experiments and section 6 gives the proof of classifi- 
cation theorem. 
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In the last section (section 7) we introduce the notion of quasi-triviality and we prove that all 
deformations are quasi-trivial. 



2 Poisson geometry 
2.1 Poisson bracket 

Definition 1 Let M be a smooth manifold, P l ^{x) a bivector (i.e. a skew- symmetric contravariant 
tensor field of type (2,0)), and f and g two smooth functions. The expression 

is a Poisson bracket on M if it defines a structure of Lie algebra on the ring of smooth functions 
on M. 

Jacoby identity is the only property of a Lie algebra that is not a consequence of the skew-symmetry 
of P lj and of the definition of Poisson bracket. 

It is well known that the Jacobi identity holds if and only if the tensor 

fipij dPi k flpki 

J l3k = {&> *h + cycUc = —P° k + — P* + — (2.2) 



is indentically equal to 0. In fact 

' n \ hi 4- rurlir = J^ h 

dx l dxi dx k 



{{L9},h} + cyclic=J^^^- (2.3) 



for any /, g and h £ C°°(M) 

Definition 2 A vector field VH associated to a smooth function H by the formula 

fjlT 

(VHy = P^— = {x\H} (2.4) 
is called Hamiltonian vector field. 

Definition 3 A smooth function f is called Casimir if {/,.} = 0. 
2.2 Schouten bracket and Poisson cohomology 

Let T l (M) be the space of i- vectors. It is well-known (see [16]) that there is an unique well defined 
R-bilinear extension of the Lie- derivative to an operator 

[., .] : T P (M) x T q (M) -> T p+q {M) 

such that 

v 

[X, A ... A X p , Q] = ^(-l) m Xx A ... A X % A ... A X p A [X,, Q] (2.5) 
i=i 
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where X k G T\M) for k = 1, .., j>, Q G H(M) and [X;, Q] = L Xi Q. 

This bilinear map is called Schouten bracket. It has the following properties: 

[p,Qi = (-ir[Q,p] (2.6) 

[P, Q A P] = [P, Q] A P + (-l)^Q A [P, P] (2.7) 

(_1)P(-D [ P) [Q ) iJ]] + (_1)5(P"1) [Q, [jR) p]] + (_l)r(5-l) [jR) [ P> Q]] = o (2.8) 

where P G T P (M), Q G r«(M) and P G T r (M). 
The last property is called graded Jacobi identity. 
It can be proved by using (2.5) and (2.6). 

In coordinates the Schouten bracket of a p-vector P and a g-vector Q is given by the formula (see 
[16]) 



[P,Q] 



k\...k 



p+q—1 



( — l) p u i f)pii-ip 1 f)nh---jg 
v L ! x k L k p+q-iQuj2---3q_ I x r fc i fep+9-i pm 2 ...i p 



p!(g-l)! »i-^2...j, 5x « (p-l)!g! * 2 -^ 1 -^ 0x u 

where 5 " is the Kronecker multi-index. 
Example 1 IfP,Qe Y 2 (M) then 



[p>QF =%r psk + ^rQ sk + (2-9) 



When Q = P we obtain 



\[P,P] = J ijk (2.10) 
and then we have the following well-known 

Theorem 1 The Jacobi identity holds if and only if [P, P] = 0. 

The last step before introducing the Poisson cohomolgy is the following 

Theorem 2 Let P be a Poisson bivector on M, then the operator dp : T q (M) — > T g+1 (M) defined 
by the formula 

d P Q:=[P,Q] (2.11) 

is a cohomology operator, i.e. d 2 p = 
Proof 

The graded Jacobi identity and the property (2.6) imply 

= [[P, Q],P] + [[P, P],Q] + [[Q, P],P] = 2 [[P, Q],P] = -2d 2 P Q (2.12) 
The last theorem allows us, following [13], to introduce the complex 

o -> r°(M) -> r x (M) -> r 2 (M) -> r 3 (M)... 

and to define the Poisson cohomology as HP*(M,P) = ker(d P ) / Im(d P ). 
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Example 2 HP 

If f is a smooth function then d P f = P %s ^, that is ifP°=Casimirs. 
Example 3 HP 1 

The cocycles are the infinitesimal symmetries (dpX = LxP), the coboundaries are hamiltonian 
vector fields. Then /7P 1 =Infinitesimal symmetries/Hamiltonian vector fields. 

Example 4 HP 2 

Let Q be a bivector. Q is a cocycle if and only if [P,Q] = 0. Therefore P + eQ satisfies Jacoby 
identity mod(0(e 2 )) . This means that the cocycles are infinitesimal deformations of the Poisson 
bracket. The coboundaries are infinitesimal deformations obtained by a change of coordinates. In 
fact Q is a coboundary if and only if Q = L X P where X is a vector field. 
Summarizing : 

iPP 2 =Infinitesimal deformations of P/Deformations obtained by a change of coordi- 
nates 

2.3 Bihamiltonian structure on M 

Definition 4 A bihamiltonian structure on M is a pair (Pi,P 2 ) of Poisson bivectors such that 
for any Ai, A2 € R, the linear combination 

A]Pi + A2P2 

is again a Poisson bivector 

In terms of Schouten bracket the bivectors (Pi, P2) are a bihamiltonian structure if and only if 

[P 1 ,P 1 ] = [P 2 ,P 2 ] = [P 1 ,P 2 ] = (2.13) 

2.4 Poisson brackets on formal loop space 

Now we want to extend the previous definitions to the loop space C = {u : S* 1 — > M}. Let A be 
the space of differential polynomials in u, that is 

/ G A ^ f = ( 2 - 14 ) 

where u (s ^ := f^. 

We observe that / is not necessarily a polynomial in u 1 . 

The role of the functions on C is played by the local functionals 

/ = J f(u(x),u x , u xx , ...)dx 

where / e A. 

1 In general / could depend on x. In this paper we will deal only with differential polynomials where x doesn't 
appear explicitely. 
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Definition 5 A (non local) multivector X is a formal infinite sum of the type 

d d 
x = x Sl '-' Sk ( y x 1 ,...,x k ;u(x 1 ),...,u(x k ),u x (x 1 ),...) - A ... A 



du Sl (xi) du Sk (x k ) 

where the coefficients satisfy the skew-symmetry condition with respect to simultaneous permuta- 
tions 

Sp, Xp < ^ 5gr, Xq 

The wedge product of a A;-vector X by a /-vector Y is defined as 
(X A F) Sl '-' Sfc + ( (xi, x k+ u u(xi), w(x fc+i ), ...) = 

_L ^ ''(*)( Ij(1)l ... )IffWl ...)yvHi) MHo( Iff(w)) ... )IffWi ...). 

Definition 6 ^4 k-vector is called translation invariant if 

X Sl '-' Sk (x u ...,x k ; u(x k ), ...) = d s x \...d s x k k X(x u x k ; ufa), ...,u(x k ), ■■■) 

where X(...) means X°~°(...) and for any t 

X(xi +t, ...,x k + t;u(xi), ...,u(x k ), ...) = X(xi, ...,x k \u{x{), ...,u(x k ), ...) 

It follows from this definition that a translation invariant multi-vector field is completely charac- 
terized by the "components" 

x *i-*k .- X (x u ...,x k ;u(x!), ...,u(x k ), ...) 
Definition 7 A k form u is a finite sum 



w — jj UJ s 1 ...s k 8u Sl A ... A 5u Sk (2.15) 



where uJ Sl ...s k 



In order to define a Poisson structure we need to introduce a criterion for the Jacobi identity. 
We have seen that in the finite-dimensional case the Jacobi identity can be written in terms of 
the Schouten bracket. 

Therefore if we will be able to define an infinite-dimensional version of the Schouten bracket we 
will be also able to define an infinite-dimensional version of the Poisson bracket. 
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2.4.1 Schouten bracket of translation-invariant multivectors 



In the case of translation invariant multivectors one obtains the formula for the Schouten bracket 
simply by translating the formula (2.9) in the new context. Heuristically this can be done by 
substituting sums for integrals, partial derivatives for variational derivatives, etc. The result is: 

Definition 8 Schouten bracket of a translation invariant p-vector P Xi -- x p and a translation in- 
variant q-vector Q Xl - Xq . 



[P, Qjzi-zp+s-i 

fepih-D'W ) { 8u(-H Ml) ) ) + 

iZ^ i , x f dQ x o(p)- x *( P+q -i)\ \ 

V \d S p x ^P+i) x aW- x a(p-l)) I _L (2 iq) 



i=0 

In the case p=2: 



^ Q]X1...X,+1 = 

'^P%(1)%(2) 



rS5„ + l S=0 \Z=1 



+ 



iZ^ 1 , N / Qp x <r(2)-X<r( q + 1)\ 

_ ( 5 s p^( 2+1 )^(i) ) _ f2 17) 



i=0 



ql \ ~°^> / V du i - s ){x a(2+i ) 



Example 5 Schouten bracket of a bivector P with the components P xy and a local functional 
I = J s i f(u(x),u x , u xx , ...)dx 

lPj] ' = L^k p - dy (218) 

Example 6 Lie derivative of a translation invariant bivector P with the components P xy along a 
translation invariant vector field Q with the components Q x 



IP, QP = 
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Example 7 Schouten bracket of two translation invariant bivectors P and Q 

1 / 8P xy 80 xy 8P xy 80 xy 

[p, Q\*y* = y( -^^d s x Q xz + -^—d s x P xz + ——d s o yz + -^—d s v P yz 

L 1 2^\8u s (x) x 8u s (x) x du s (y) y du s (y) v 

apzx 80 zx apzx 80 zx 8P yz 

+_ 8 s O zy + — - 8 s P zy + — 8 s O xy + — - 8 s P xy + — 8 s Q yx + 

+ du*(z) M + du*(z) z + du*(x) ^ + du'(x) x + du°{y) yUi + 

8Q yz 8P yz 80 yz \ 

Remark 1 

The operator is the usual partial derivative when it acts on functions depending on x while 
it is, by definition, equal to when it acts on functions not depending on x. In other words 

du{x) 



du(x) 
du{y) 
du(x) 



= 



Remark 2 

The formula (2.19) can be obtained by the formula (2.17) for q — 1 by using skew-symmetry of 
P. 



Definition 9 A translation invariant Poisson bivector P is a translation invariant bivector such 
that 

[P, P] Schouten = 

As in finite dimensional case a translation-invariant Poisson bivector P xy defines a Poisson struc- 
ture on the loop space C. The Poisson bracket of two local functionals p ,/ 2 is given by the 
formula 

J s iJ s iSu(x) P 5u(y) dxdy 

To define the Poisson cohomology on the space of translation invariant multivectors we have to 
prove that the operator dp := [P, .] associated to a Poisson bivector P satisfies the condition 
dp — 0. We have seen that this condition is satisfied (for a Poisson bivector) if the graded Jacobi 
identity holds and moreover that the graded Jacobi identity follows from (2.5) and (2.6). 
The property (2.6) is obvious. As far as it concerns (2.5) it will be sufficient to analyse a particular 
case to understand how the proof works in general. 

We want to check the "Leibniz rule" (2.5) when p = 2 and Q is a bivector. In this case we have 
to prove: 



(2.21) 



[X A Y, Q] = Y A [X, Q] - X A [Y, Q] 



(2.22) 



Left hand side 

By using formula for the Schouten bracket of two bivector, we obtain 



[X A Y, Q] 



xyz 



(dX(x 



\du s (x 
(dY{x 



\du s (x 
dX(z 



du s (z 
dY(x 



du s (x 
fdX(y 



\du s (y 
dY(z 



du s (z 



Y{y)-X{y) 
X(y) - Y{y) 
Y{x) - X{x) 
X(z) - Y(z) 
Y{z)-X{z) 
X{y) - Y(y) 



dY(x 



du s (x 
dX(x 



du s (x 
dY(z 



du s (z 
dX(x 



du s (x 
dY(y 



du s (y 
dX(z 



du s (z 



(d s x P xz ) + 
(d s y P yz ) + 

(9 S z P Z y) + 
(d S x P X y) + 

(d s z P zx ) + 



dR 



xy 



du s (x) 



dR 



xy 



du s (x) 

0P ZX ' 
du s (z) J 

dP zx 
du s (x) J 

OP 

ul yz 

du s (y) 



dR 



yz 



du s (z) 



mX)Y(z)-X(z) (d s x X)) 
{(d s y X)Y(z)-X(z) {d s y X)) 
)mX)Y(y)-X(y) ftX)) 
)mX)Y(y)-X(y) (d s x X)) 
{(d s y X)Y(x)-X(x) (d;x)) 
mX)Y(x)-X(x) (d s z X)) 



This formula can be written collecting terms in Y(x), Y(y), etc.: 



Y(x) 



dX(y) 



du s (y) 
Right hand side 



du s (z) 



du s (y) 



du s {z) 



(2.23) 



Y A [X, Q] - X A [Y, Q] = 

= \ E(- 1 )* an(ff)y (^i))[^ p ]w-(.) - \ E(- 1 ) S9ri(CT)x (^))[ F ' p ]^( 2 )^(3) = 

= l - (Y(x) [X, P] yz - Y(x) [X, P] zy - Y(y) [X, P] xz + Y(y) [X, P] zx + Y(z) [X, P] xy 
-Y{z)[X,P] yx )- l -{Y~X) 
Let us consider the first two terms: 

l -{Y{x)[X,P] yz -Y{x)[X,P] zy ) 

They can be written as 
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By comparing this expression with (2.23) we realize that they coincide because of the skew- 
symmetry: 

dP zy _ 0P yz ^ 2 24 ^ 



du s (y) du s (y) 

dP zy _ dP yz (2 25) 



du s (z) du s (z) 
Analogously for the other terms of the left and right hand side. 

2.4.2 Local multivectors and their Poisson cohomology 

Definition 10 A local k-vector is a traslation invariant k-vector such that its dependence on 
Xi,...,x k is given by a finite order distribution with the support on the diagonal X\ = ... = x k . 

In coordinates a multivector X has the form 

X = X(u(x 1 ),u x (x 1 ),...)5^\x 1 -x 2 )...5^\x 1 -x k ) (2.26) 

P2,P3,-,Pk>0 

It is easy to check that P xy = u5'(x — y) + \u x 5(x — y) is a local bivector. In fact 



D yx = u(y)5'(y -x) + ^u y 5(y - x) = -u x 5'(x - y) - u x 5(x - y) + ]-u x 5(x - y) = -P xy 



But is not equal to -|^r = 0! 

au(x) ^ ou(x) 

This problem can be solved by writing P xy in a suitable form. If we write P xy as P' xy = P xy = 
~2^Pxy ^yx) then 

Qpl Qpl 

xy - - yx (2.27) 



du{x) du{x) 

This is true in general: if we want to use the same formulae valid for non-local multivectors we 



,nuie wini wits upeiaiuis ^ ^ 

rule is to write the multivector T xl '-> Xn in the form 



have to write the local multivector in a form "compatible" with the operators du s^ ■ The practical 



J_ \^ fiygn(a-)rpx^ {1) ,...,x a(n) 

n\ ^ 

Analogously to the non-local case the formula (2.17) allows us to define a cohomology operator 
dp starting from a local Poisson bivector P: 

n r dp r l dp r 2 d P F 3 dp k , , 

U ~ 1 ~ 1 ~ 1 local ~^ 1 local —>■■■—> 1 local--- {Z.ZO) 

where T° is the space of local functionals, T 1 is the space of non local translation-invariant vector 
fields and Tf ocal is the space of local translation-invariant A;- vector fields. It is easy to check that the 
cohomology groups we have introduced in this way have the same meaning that the cohomology 
groups of finite-dimensional Poisson manifolds (see the examples in the section devoted to Poisson 
cohomology of finite dimensional Poisson manifolds). 



10 



3 Deformations of bihamiltonian systems of hydrodynamic 
type 

In this paper we deal with Poisson bivectors of the form 



P*y = <f>( u )5W(x - y ) + i^M* - y) + e " p S 



(3.1) 



k 



where 



k+l 



(3.2) 



s=0 



where A^ s are differential polynomials. Introducing the following gradation in the space A of 
differential polynomials: 



The part of order is called local Poisson bracket of hydrodynamic type (see [6]). 
Analogously to the finite dimensional case we have the following 

Definition 11 A bihamiltonian structure on the loop space L is a pair (Pi,P 2 ) of local Poisson 
bivectors such that for any Ai, A2 G R, the linear combination 



is again a Poisson bivector. 

Definition 12 An evolutionary equation 

u t = F(u(x),u x ,u xx , ....) 
is called Hamiltonian if it can be written in the form 

u t = {u,H} 

where H is some local functional. 

Definition 13 An evolutionary equation is called bihamiltonian if and only if it is hamiltonian 
with respect both Poisson bivectors of a bihamiltonian structure. 



deg(f(u)) = 
deg(u {k) ) = k 



we require also that 




(3.3) 



A1-P1 + A2-P2 
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Example 8 Reseating x — > ex the KdV equation u t = uu x + u xxx one obtains 

u t = e(uu x + e 2 u xxx ) 

One usually introduces slow time variable t — > et to rewrite the last equation into the form 

u t = uu x + e 2 u xxx 

This is small dispersion expansion of the KdV equation. It is a bihamiltonian equation. 
The first Poisson bivector ( Gardner- Zakharov-Faddeev bivector) 

P*y = S'(x - y) (3.4) 
and the second Poisson bivector (Magri bivector) 

{u(x),u(y)} 2 = u5'(x -y) + ^u x 5(x - y) - e 2 5 3 (x - y) (3.5) 
belong to the class of bivectors introduced before. 

It is easy to check that the pair (P 1 ,P 2 ) of the Gardner-Zakharov-Faddeev bracket and of the 
Magri bracket is a bihamiltonian structure. It is called Magri bihamiltonian structure. 

Definition 14 The group of transformation 

u -> u = ^2e k F k (u,u x ,...) (3.6) 

k 

where F k e A, deg(F k ) = k and ^ ^ is called Miura group. 

Degiovanni, Magri and Sciacca (see [4]) and Getzler (see [11]) solved independently the problem 
of studying the action of the Miura group on the bracket (3.2). 

More precisely they study the following problem: does it exits an element of the Miura group that 
transforms the bracket (3.2) into the bracket (3.4)? 

This problem can be reduced to a cohomological problem. In fact we have seen that it corresponds 
to the study of the second group of Poisson cohomology associated to the bracket (3.4) 2 . 
HP 2 (P, n)=Infinitesimal deformations of P/Deformations that can be obtained by 
infinitesimal change of coordinates of the form (3.6). 

This cohomology group is trivial (see [4] and [11]). 

In this paper we deal with deformations of bihamiltonian structure of hydrodynamic type. 

Definition 15 A deformation of order k of a bihamiltonian structure of hydrodynamic type is a 
pair (Pi,P 2 ) of the form (3.2) such that Pi — \P 2 satisfies Jacoby identity for every A up to the 
order k: 

[Px, Px] = [P 2 , P 2 ] = [Px, P 2 ] = o(e fe ) (3.7) 

From the triviality of second Poisson cohomology group it follows that we can always assume one 
of the brackets of the form (3.4). 

In the next section we explain how to classify these deformations. 

2 in our case the role of infinitesimal change of coordinates is played by the Miura group 
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3.1 Classification of deformations of bihamiltonian structure of hy- 
drodynamic type 

We are interested in the following question: which deformations of a bihamiltonian structure 
of hydrodynamic type are trivial? In other words which deformations can be obtained from a 
bihamiltonian structure of hydrodynamic type by the action of Miura group? 
We start considering first order defomations, that is Pi = = 5 1 (x — y) and P 2 = P^ + eP 2 ^ 
By definition we have 3 

[P 1 ,P 2 ] = o{e) (3.8) 
[P 2 ,P2] = o(e) (3.9) 

The equation [Pi, P 2 \ = o(e) implies [Pf > \ P^] = diP 2 ^ = and from the triviality of the second 
Poisson cohomology group it follows that there exists a vector X 2 ^ such that 



P« = d ± X^ (3.10) 

The equation [P 2 , P 2 ] = o(e) implies that 



[? 2 (0) , P 2 (1) ] = d 2 P 2 {1) = d^X^ = -d^X™ = (3.11) 
where we have used the fact that (di + d 2 ) 2 = 0. 

Among all deformations that satisfy the equation d\d 2 X 2 ^ = we have to select trivial deforma- 
tions, that is the deformations that can be obtained by infinitesimal change of coordinates. In our 
case this means that there exists a vector field X such that 



LiexPr = (3.12) 
L^P 2 (0) = P 2 {1) (3.13) 

Theorem 3 A deformation P\ = P^ — X(P 2 °^ + ediX 2 ^) is trivial X 2 ^ = d±a + d 2 b where a, 
b are local junctionals. 

Proof 

= d ia + d 2 b P 2 (1) = d x d 2 b = -d 2 d x b 
This means P^p = Lie^pf^ with X = —d\b. 
Moreover Lie^P^ = 0. 

Lie^Pf^ = =>■ X = d\b (the first Poisson cohomology group is trivial). 
-d x d 2 b = d 2 d x b = LiexP 2 0) = P 2 ] = d^xf* => x ? = -d 2 b + d x a 



3 It is well-known that [Pi, Pi] = 
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From the last theorem it follows that the elements of the group Ker{d\d<2)l (Im(di) + Imfa)) 

are the non trivial first order deformations. 

Remark 

For higher order deformations we can repeat the same arguments and obtain the same result. 
Consequently, in general: 

1) A k order deformation P^ can be represented in the form 

pf _ AP (o) + iPf = P 2 (0) - APf + ^ X 2 ] (3-14) 
i=i i=i 

where due to our definition of a gradation in the space A of differential polynomials (see (3.3)) 
we have necessarily deglX^) = i (see the form of the formula (2.19) when Pi — 5'(x — y)). 

2) P 2 (fc) is trivial if and only if xf ) = d x A + d 2 B. 

Remark 1 Also in this case "trivial" means that P^ can be eliminated by the action of Miura 
group. Obviously the triviality of P^ doesn't imply that the /c-order deformation is trivial. 

Remark 2 

From the (2.18) it follows immediately that d± and <i 2 increases the degree of a local functional by 
one 4 . Then the degree of A and B must be equal to k — 1. 

In the next section we summarize the results of the classification of deformations up to fourth 
order. 

3.2 Classification of deformations: results 

Theorem 4 Up to the fourth order all deformations of a bihamiltonian structure of hydrodynamic 
type (see definition 15) can be reduced, by the action of Miura group to the following form: 

u5 {1 \x -y) + ^uj(x -y)- \5 {l \x - y) + 

+e 2 (~2s5 3 (x -y)- 3d x s5 2 (x - y) - d 2 x s5\x - y)) + 

+e 4 (-2S5 5 {x - y) - 5 (d x s) 5\x - y) - 10 (d 2 j) 5 3 {x - y) - 10 (d 3 j) 5 2 {x - y)+ 

-3 (<9 4 §) 5\x -y) + 2w5 3 (x -y) + 3 (d x w) 5 2 (x - y) + (<9» 5\x - y)) + 0(e 5 ) (3.15) 

whith 

ds 

W = 2 g^ U xx (3.16) 

where s is an arbitrary function of u and s = — 2s|^. 

From this theorem it follows immediately: 

4 we will identify the degree of a local functional f fdx with the degree of its density / 
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Corollary 1 Up to the fourth order every deformation of Magri bihamiltonian structure is trivial 
Proof 

Indeed, in this case s = constant ^ and then s = 0, w = 0. 



4 The deformed hierarchy 
4.1 Integrals of motions 

As we have already been said bihamiltonian structures give rise to an infinite number of " almost- 
constants" of motions. From these constants one can construct an hierarchy of hamiltonian equa- 
tions. 

In this section we show explicitely how to produce this hierarchy and we check that the corre- 
sponding flows commute up to the order of the deformation. 

The technique is well known for the KdV equation (see [9]): one looks for a solution of the equation 

Pxv = (4.1) 

in terms of a formal series: the coefficients of 1-form v are variational derivatives of the integrals 
of motions. 5 

We apply the same technique to the second order deformations (see classification theorem). 

We start with the equation P\v = 0. This equation implies vP\v = 0. From skew-symmetry of 

Pa it follows that = vP\v = d x (...). That is (...) = /(A). More precisely 

v(-Xd x + ud x + \i x + e 2 (2s(u)dl + 3(d x s)d 2 x + (d 2 x s)d x ))v = (4.2) 
From (4.2) it follows: 

9x{^v 2 (u - A) + e 2 (2svv xx + s x vv x - sv 2 x )) = (4.3) 

that is 

l -v 2 [u - A) + e 2 {2svv xx + s x vv x - sv 2 x ) = /(A) (4.4) 
By choosing /(A) = — | we look for a solution of the form 

i=0 

where obviously po — 1 . 

By straightforward calculation we get 

1 PiPi f U 2/r. \ 1 m 

2 2^p = 2^ (2 PiPj + € ( PiPjxx + SxPiPjx ~ SPixPjx >\^ ( ' 

i+j=k i+j=k—l 

5 this is equivalent to the exactness of the 1-form v with respect to the vertical differential of the variational 
bicomplex (see [3]). 
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for k — 1 (4.6) implies 



that is 



for k=2: 



^(PoPi + PiPo) = ^upl (4.7) 
Pi = \ (4-8) 

P2 = ^u 2 + ^{2su xx + s x u x ) (4.9) 

and so on. We observe that the coefficients will have always the form 

Pl = f(u) + 0(e 2 )+0(e 4 ) + ... (4.10) 

where f(u) is a polynomial in u. 

Now we want to prove that the coefficients pi are, up to the second order variational derivatives. 
We consider a curve u(t) on the loop space L. By differentiating the equation 

^v 2 (u - A) + e 2 (2svv xx + s x vv x - sv 2 x ) + ^ = (4.11) 
along the vector field it, tangent to this curve, we get 

(v v 
4:sv xx + 4s-v xx + 4:sv' xx + 2s' x v x + 2s x -v x + 2s x v x + 
V ' V 

V V 

By straightforward calculation one can write the last equation in the form 

iiv = d t (-2^) + e 2 (d t (-4sv xx - 2s x v x ) + d x [^s\^j + 2^ (sv x - s x vfj (4.12) 



This equation implies 



J u V dx = j t (^j fdx^j+0(e 4 ) (4.13) 



for some / = f(u, u x , u xx , ...). 

In fact from (4.10) it follows that v has the form 

v = v ° + 0(e 2 ) (4.14) 

where v is a function of u and 



sv° x - s x v° = 



Now we can formulate the 
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Theorem 5 p { = g + 0(e 4 ) /or some junctionals Ii . Moreover 

{IiJih = 0(e 4 ) (4.15) 

Proof 



l* vdx =W, fdx ) +o{e4) =LX 



*4r + ••• ) dx + 0(e ) = ... 

ou OU T OU TT 



By integrating by parts, we get 



9f a fdf\ , a2 f df \ , 



Then 



t. = f^ + 0(e 4 ) (4.16) 

OU 



that is 

5/, 



Pl = ^ + 0(e 4 ) (4.17) 



for some functional /«. 

Moreover, by definition, the coefficients Pi of 1-form v satisfy the equation 



(p 2 - apo (^E |) = ( 4 - 18 ) 



In terms of L this condition can be written as 



, 5m 7 A* ^ V 5u J X 

t=0 x 7 i=-l v 7 



From this equation it follow immediately: 



5u 


= 


(e 4 ) 


p 6I i 

2 5u 


= Pi 


5u 



By using these Lenard recursion relations it is easy to prove the therem (see for details [9]) 
4.2 Soliton solutions 

Now we study one of the equation of the hierarchy. More precisely we concentrate on the equation 

5I 2 3 2 / <9s 1 d 2 s 3 



u t = d x — = -uu x + e ( s(u)u xxa! + 2—u x u xx + ^-Q^Ux ) ( 4 -20) 



When s = constant this is KdV equation. 
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4.2.1 One soliton solutions 

We look for solutions of (4.20) of the form u(x, t) = u(x + ct) = u(z). 
By substituting we get 

3 2 / ds 1 d 2 s o\ . . 

cu z = -uu z + e \s{u)u zzz + 2—u z u zz + -—u z \ (4.21) 

This equation can be written as 

d z (cu - ix 2 - e 2 (d z {su z ) - l -(d z s)u)^ = (4.22) 

that is 



n, - ~ir -c 2 [d z (su z )-^(d z s)u z ) r-j (1.23) 



By multiplying for u z we get again a total derivative 



that is 



0_[^cu 2 -^-e 2 (\su 2 z )-c x u)=() (1.21) 



1 1 / 1 

-cm 2 - -m 3 - e 2 f -sm 2 ) - ciu = c 2 (4.25) 



This equation can be written as 

,du 



where 



(=) 2 = F(u) (4.26) 



To obtain the solution one has to invert the following integral 



z - z = ± 



1 Tw*** (4 ' 28) 



4.2.2 Case s(u)=u 

In this case F(u) = where P(u) is a polynomial of degree 3. It is well known (see [5]) that 
one soliton solutions occur when F(u) has one simple zero and one double zero. 
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By using the formula (see [15]) 

x -a 3 



(a — x)(x — a 1 )(x - a 2 ) 
where a > u > ai > a 2 > a 3 and 



dx 



2(a - a 3 ) 



-n 



((a -a 2 )( ai -a 3 ))2 V a x - a 3 



arcsm 



[a x - a 3 )(a - x) \ 2 



A: 



(a - ai)(x - a 3 ) 
a - ai)(a 2 - a 3 



ai - a 



(a - a 2 )(ai - a 3 ) 

n(0,i/,fc) = r 

Jo (1 — vsin 2 (j)) (1 — k 2 sin 2 (p) 2 
we obtain, when a\ = a 2 and a 3 = (see [1]) 



(4.29) 



(4.30) 

(4.31) 
(4.32) 



-2 - Zq 



±- 



2ai 



/ 



(ai (a - ai)) 2 



-2/n 



V 



^ / ai(a --uA 2 , 
\^u(ao-ai) y 

T 

/ ao(M-ai) \ 2 
\ y«(ao-ai) y / 



(4.33) 



The speed c of the wave and the constants of integration ci and c 2 can be easily written in terms 
of the coefficients Qq, a±, a 2 , a 3 : 



c = -(a + 2ai) 
ci = i(2a ai + a?) 



c 2 = --aoOi 



(4.35) 
(4.36) 
(4.37) 



The expression (4.34) can be inverted numerically. To test the existence of two solitons solutions 
we have used two such solutions with different speed as initial condition. The results of numerical 
experiments are described in the next section. 
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5 Numerical experiments 



0.0019 



0.0012 




20 40 60 80 100 120 140 160 180 

In this section we analyse the equation (4.20) for s(u) = u, i.e. 




We write this equation in the form 

ut + F x = 

where F = — j 2 -. To make numerical experiments we have used a two-steps Lax-Wendroff scheme 
(see [14]). This scheme is characterized by an auxiliary step of calculation 

where F? +1 = F(u] +1 ). 
The main step is 

u n + l = - ^{F n +;> - F n+ h 

The first figure illustrates two solitons of different height and speed we have used in numerical 
experiments. 

We have choosen periodic boundary conditions. 

The results of our experiments is showed in figure 2. The waves move from the right to the left. 
The three pictures of the second figure show the solitons before, during and after the collision. 
Like the soliton collision in integrable systems they reemerge with the same shape. 
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0.0017 - 
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0.0011 - 
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Remark 

Due to numerical instability for big amplitude we performed numerical experiments with small 
and slow waves. 

We don't know if the origin of this numerical instability is the alghoritm we have used. 

6 Classification of deformations: proof 
6.1 First order deformations 

We have seen that in this case the non trivial deformations are the elements of the group 
Ker{d x d 2 )l (im(di) + Im(d 2 )). 

We start looking for the solutions of the equation did 2 X = 0. 
First of all we calculate d 2 X 




+ 



-u y 5(y-x)) 



du s (y) 



dX 



where, by definition, 5'{x — y) — 5^(x — y) — d x 5(x — y). 
Moreover by using the formula 




(6.1) 
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we get 



and 



\ {d^{x) {u6 ' {x ~ y) + \ uJ{x ~ y)mx) ~ d^y) {udAv ~ x) + \ uAv " x)){d y x) 

= l - (xS'(x -y) + \d x X5(x -y)- Xd y S(x - y) - \d y X5{y - xfj = 
= X5'(x -y) + ^d x X8(x - y) 



dy(ud y S(y -x) + -u y 5(y - x)) ^ - d^(u5 (1 \x - y) + -uj(x - y)) 



1 f)V 

= <? y (-u(x)5M(x -y)- -uj(x - y))^ + 

1 dX 

-d s x (-u(y)d y 5(y - x) - -u y 5(y - x))- 



(y) xx v »> 2 x v *"0u fl (x) 
dX 



2 y v * yy <9w s (x) 
(-1Y +I (u(x)5^ +I \x -y) + -uj^(x - y))£^ + 



-u{y)5^ +l \x-y)+ l -u y 5^\x-y)y 



OX 



2 y v ^du^x) 
Summarizing we obtain 



-, s+l 

{d 2 x) xy = xsW(x - y ) + -d x X5(x -y) + J2J2 c ^ {s+1 ~ q \x - y) (6.2) 



s q=0 

with 



1 ( S 
+ - 



2V 9 -iA- 1),+1 "4 9 *i^)) + a^)" W ) < 6 ' 3 > 

In this case deg(X) = 1 i.e. X = s{u)u x . This implies that we can write the sum (6.2) as 

(d 2 X) xy = 

-d x X + c 2 i + cio <J(x - y) + (X + c n + c 00 ) (x - y) + c i5 (2) (x - y) (6.4) 
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with 



d 2 s 1 ds\ 2 ( ds 1 
<9 2 s 9s\ 2 ds 



en = 2u—u x — su x 



ds 

du 
ds 

c 00 = -2u—u x 
du 

c i = 

Substituting the last equations in (6.4) we obtain (d2X) xy = 0. 

This means that all fields X of degree 1 belong to ker[d\d2). 

Now we have to calculate the trivial field i.e. the fields X = d\A + d 2 B. 

where 



is 1 

5 = / B(u)dx 
Js 1 



Using the formulae 



and 



we get 



X A 

diA = -d x — (6.5) 



d 2 B = -d x [ u 5 ^ + \^u x (6.6) 



. , . . d 2 A(u) IdBiu) d 2 B(u) . 



This shows immediately that all deformations are trivial. 



6.1.1 Explicit form of the deformations 

We have to calculate d±X for an arbitary field of degree 1. 

d 2 X = implies d\d 2 X = —d 2 d\X = 0. Then there exists a vector field X of degree 1 such that 
d x X = d 2 X = 0. 

This argument can be also used as an alternative proof of triviality of first order deformations. 
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6.2 Second order 

Now we consider the deformation P 2 = P 2 ^ — XP^ + € 2 P 2 2 ^ 6 . 

Also in this case the non trivial deformations are elements of the group Ker{d 1 d 2 ) / (Im(di) + Im(d 2 )). 
We start again considering the solutions of the equation d x d 2 X = 0. 
In this case the sum (6.2) becomes 

(d 2 X) xy = (^d x X + cio + c 2 i + c 32 ^ S(x - y) + (X + c 00 + c u + c 22 ) - y) + 

+ (c i + ci 2 ) 5 (2) (^ - y) + c 02 5 {3) (x - y) 
where X = squ xx + S\U 2 X and 

cio = -a a I 

C21 = ^ W x ) ~ 2d^ U ™ + 2 dx W x ) Ux 
a3 f dX \ IdX 1 a2 ( dX } 

C32 = ~ udx ) ~ 2d^ Uxxx ' 2 d * 1 u 

« 9X 

c 00 = -2u—- 

ou 

. , dX\ Ox 

Cn = 2lXd x - 7^ M x 



a; 



du x J du x 



C22 = ~ 3ud * ( 15) " ~ 9x (J^r ] u - 



, dX \ dX dX 
coi = 0, C12 = -3ud x - - 3- u x , c 02 = -2u 



f) J f) x " UZl £) 

I XX XX 



Hence we can write 



(d 2 X) xy = J2hS (k) (x-y) (6.8) 



fc=0 



with deg{b k ) = 3 — k. By straightforward calculation we get 

6 The trivilaity of first order deformations means that we can always kill the term of first order in e with a change 
of coordinates 



24 



b = aU xxx + bu x U xx + Cu\ = ( ~ 2u ^ + 2uS l J U xxx + 

s - iu d2 ' So + 4m— - — ^ u u + (u d2 ' Sl + -— - u——— - u 3 
du 2 du du ) \ du 2 2 du du 3 2 du 2 J x 

h = du xx + gu 2 x = (^~ 5m ^ + 4ms i - s o^ u xx + (^ u ^ ~ 3m ~^ ~~^~ Sl ^ u2 

02 = = — O M— h Sq Ma 



du 



I = — 2uSr 



Now we can calculate explicitly the equations did 2 X = in terms of the coefficients 6 , 6i, 6 2 , 63. 
By using the formula 



(did 2 ) Q 



1 y d(d 2 X) ) (t+1) _ _ 1 ^ 0(rf 2 X) ) (m) _ 



2^ ^(x) ' 2 U dut (y) 



I d(d 2 X) z (t+1) _ g(^ r( W 
2^ du*(z) 2^ 

1 A WOp ^), _ x) _ 1 f gWOg ^Df 



2; — x 



that is 



4 g (all) ) * - »>* m, <* - *> + ^ (3, <* * - * + 



In order to obtain a sum where appear only terms with S^(x — y)5^\x — z) and where the 
coefficients depend only on x we use the identity 



5(x — y)S(x — z) = 5(y — x)5(y — z) = S(z — x)5(z — y) (6.9) 
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and formula (6.1). For example we can write 

4£(jy*-^ (,+i) ^)= 

In this way, after a long but elementary calculation we obtain 

d x d 2 X = 2f5 A 5 3 + A 5 2 + 3/ x <5 A 5 2 + A 5 1 (6.10) 
where A 5 j = \ (8®(x - y)5^(x - z) - 5^(x - y)8®(x - z)) and 

/ = ^r- 9 *k- -^r + ^hr ( 6 - n ) 



ou x \ ou xx J ou \ ou 

Therefore the equation d^X = is equivalent to the equation / = 0. By substituting &i, 6 2 , 6 3 
in this equation we get the condition 



6.2.1 Trivial deformations 

The differentials di and <i 2 increase the degree by one. Because the degree of X is 2 the trivial 
deformations can be written as X = d\A + c?2-B = where the degree of the densities of A and B 
is 1. But the variational derivative of local functionals with densities of degree 1 vanishes: 

s U m^ _ ^ { _ m ( V\ = a, _ = 



5u ^ v ' x \ du^ Ou 

i=0 



Then all the deformations d\X\ (with X\ = s u xx + fpw 2 ,) are not trivial. 
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6.2.2 Deformations: explicit form 

By using the formula (2.19), we get 



^ 5 ' {y ~ d^) {sUyy + d^ Uy) " {dSJ ' {x ~ y)) d^) {sUxx + o^ ul) ) 
(t> + U Uy ) 6W{y - x) - (th* + U ul ) 6(1){x - y) + 

^n)j 6®(y -x)- (2^) 5^\x - y) + s(y)5^(y - x) - 8 (x)6®(x - y) 



By observing that 

ds d 2 s 2 o2 

and by using the identity (6.1) it is easy to get the formula 



= d ± xP = -2s5 3 (x -y)- 3d x s5 2 (x - y) - d 2 x s5\x - y) (6.13) 
6.3 Third order 

The condition of compatibility [Pi,P 2 ] = implies that P 2 = diXf 1 and the Jacoby identity 

[P 2 ,P 2 ] = implies that d 2 P 2 (3) = -cW 2 x[ 3) = 0. 

Remark 

There are no conditions containing the second order deformations. 
We start again calculating d 2 X. From (2.19) it follows: 

(d2X) xy = 

QdjJT + cio + c 2 i + c 32 + c 43 ^j S(x - y) + (X + c 00 + c u + c 22 + c 33 ) <5 (1) (x - y) + 
+ (coi + cia + c 23 ) 5 {2) (a; - y) + (c 02 + c 13 ) 5 (3) (^ - y) + c 03 5 (4) (a; - y) 
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where X = squ xxx + s\u x u xx + S2U X and 



Then we can write 



do = ~d x (u^j 



o2 fdX\ IdX 1 . /ax\ 
n , / ai \ 1 ax 1 ^ 9 / dX 



o,3 / «9X 



1 dX 
2du 



c 43 = 



ax 

c 00 = -Zu— 
ou 



*\du xx ) 2du xx xxx 
/ 91 \ 1 ax 1 3 / c)X ' 

I ~p. I U xxxx + —U x O x I — 

\ OUxxx J A OU xxx £ \ ou xxx _ 

r) Y 



cn = 2^ 



ax 



a-u. 



C22 = -3w<9^ 



XX 



c 33 = Audi 



x 1 0u 



XXX 



Coi = 

C12 - 3ud x 



dX 

du x 



ci 3 = 4^ 

C 3 = 



9u 



XXX , 



-I/, 

ax 



ax 



\du 

ax \ 5 ax 



2 „ ^aia: u x d x 

OU xx 



' ax 
. a« 



';rx 



2a« 



-It 



XXX 



-Mo 



.xx + ~U x dl 



' ax 

. du 



'XXX , 



C23 = 6udx {du~r2du~ 

\ UU XXX/ ^ UUj XX: 

n ( ox \ ax 



"""XXX 

ax 
'f> — Ux 

W'xxx 



Mix ~\~ ^l^xdnc 



' ax 

. du : 



XXX . 



(d 2 X) xy = J2b^ k \x-y) 

k=0 



(6.14) 
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with degipk) = 4 — k and 

b = au xxxx + bu x u xxx + c (u xx f + du 2 x u xx + g?4 = 



= + 



-3w-^ + + 6ms 2 J {u x u xxx + (w^) 2 ) + 



-to-Qj + Uu— --— + 3s 2 + 6«— + -_ ) u x u xx 4 



V " ^ 2 ou 2du 2 du 4 2du 3 U du 3 J ^ 

h = hu xxx + lu x u xx + mu x = 

/ <9s 3 \ / 9«i <9 2 s 3<9s \ 

= I 2m— - «si - -s I + I I2us 2 - 3si - 9m — + 12u ~q^2 + 2 J + 

. ( A ds 2 d 2 s 1 ds! d 3 s 3d 2 s \ 3 

&2 = P«:r:r + ?M 2 = 



<9s 9 \ / <9si <9 2 s 3 ds ( 



- 1 ^ 6U fa - 2 So ) + {- 3u fa ~ 351 + ^ + 2fa ]> ' 

h_, = /,•// , = ( -2usi + 1//^- - 3s j ", 
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6.3.1 Deformations 

In terms of coefficients bk the equation (d\d2X) xy = can be written simply by adding to the 
equation obtained in the previous case the terms 

1 ( db 



2 \ dn 



2 \dUyyyy 

1 / db 2 



2 \^ C^Uyy 

+ lf db 



5(x - y)5 {5 \x -z) + \ (-^-) 5 {1 \x - y)5 {A \x - z) 

l \ ou xxx J 

5^(x - y)5V>(x -y) + \ 8<®{x - y)5^(x - z) 

) 8(y - x)6<®(y - z) - \ {^) - ^ (A \v ~ z) 
) 5^(y - x)6®(y - z) - \ 5^(y - x)5^(y - 



2 \du 



5(z - x)5^ (z-y) + 1 - ^ (* ~ ^ (* ~ V) 

-\ iw^A ^ - z ^ x -y>-\ brA s(1) ( x ~ z ^ x ~ 

L \ OU xxxx I Z \ OU xxx I 



L xxxx 

1 / dh 



J'2 



+ 



2 V d'ti/j^Q^ 
1 f db 



2 \ 8v 

Z, \ VLlyyyy 



1 ( db 2 



2 \dUyyJ 

1 / <%o 



2 \<9lf,2ZZZ 

1 / db 2 



) ^(s - ^ 3 )(* - y) - i (J^) *< s >(* - *)<j( 2 >(* - y) 

-) S(y - z)5^\y -x) + \ (^L.) S^(y - z)5^\y - x) 
v J 1 \ au yyy / 

) *< 2 >(y - z)6®(y - x) + \ 6<*>(y - z)&\y - x) 

) ^ - y)*< 5 >(* - x) - \ (J^j 6W(z - y)6W(z - x) 
6W(z - y)6W(z -x) + 1 - 6®(z - y)SM(z - x) 



2 \du z 

By introducing the function / and g defined by the formulae 



du x \ du xx J x \ du xxx } du \ du 
db\ db 2 db 3 



we can write the result in the form 

d x d 2 X = k 6i A 6j ( 6 - 15 ) 



30 



where 



fo5 


= /23 = 2<? 


fl4 


= 5g 


/04 


= %z 


/l3 


= &9x 


/03 


= 2f + 4g xx 


fl2 


= 3/ + 


fo2 


3fx~^~ 9 XXX 


foi 


fxx 



Therefore the equation d\d 2 X = is equivalent to the equations / = and # = 0. 

By substituting the coefficients bi, b 2 and 63 in these equations we obtain that the last equation 

is identically satisfied and the first is equivalent to the condition 

2s2 -^7 + 5^ (6 - 16) 



6.3.2 Trivial deformations 

In this case trivial deformations are d\A + d 2 B with 

A= (A (u)u xx + A^ujuDdx 

B= (B (u)u xx + Bi(u)u 2 x )dx 
Js 1 

By using formula (2.18) we get 

d x A + d 2 B = (-2 - Ai \ - 2u - B 1 ) ) „,., ,. - 



| d 2 A dAA ^ f d 2 B _dB 1 \_ fdBo_ B \\ + 
1 du 2 du ) \ du 2 du ) \ du J J 
,^o_^i\_ u fd Z B d 2 Bi\ 1 fd 2 B dB 1 



du 3 du 2 J \ du 3 du 2 J 2 \ du 2 du 
If we call A:=A 1 -^ and B:=B 1 - we can write 
diA + diB = 

( 7 f A dA A dB ~\ (d 2 A A d 2 B dB 

(2A + 2uB)u xxx + 4— + 4u— + B \ u x u xx + — + 4u—r + — 
V / \ du du I \ du 1 du 1 du 

Now we can prove that all deformations P 2 are trivial. 
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Proof 

The trivial deformations satisfy equation (6.16). In fact 



d 2 A ^d 2 B 8B\ d ( .8 A 
du 2 du 2 du ) du 



(dA dB ~\ d 2 / \ 
4— + 4m— + B + — ( 2A + 2uB) = 
OU OU J OU 2 \ J 



Every field X = squ xxx + s\u x u xx + s 2 u\ with coefficients so, si, s 2 satisfying equation (6.16) can 
be written as d\A + d 2 B by choosing 

~ S U f 8sa 



2 3 V du 
6 1 fto ds ° 



and this choice is always possible. 



6.3.3 Deformations: explicit form 

By using the formula (2.19) 

( d 
d x X = (i d y 5 '(y - x )) g u (s)( y ^ s ° u yyy + s ^ u y u w + S 2 U D+ 

d \ 

- (9 s x S'(x - y)) du(s) ^ (s u xxx + s x u x u xx + s 2 ul)J 

and condition (6.16) it is easy to get the formula 

PP = dxXP = -2t5 {3) (x -y)- 3d x t5 {2) (x -y)- d 2 x t5 {1) (x - y) 
where t is an arbitrary differential polynomial of degree 1. 

6.4 Fourth order 

In this section we consider the case: P 2 = Pp + £ 2 Pp + ^P^ - 7 

The compatibility condition implies P 2 ^ = diX^ and the Jacoby identity [P 2 , P 2 ] = o(e 4 ) implies 



d\d 2 x\ — - 



d\X<^ i d\X2 



= 



7 The third order term can be always killed by changing coordinates and the second order term is 
P 2 (2) = d.X^ with X< 2 > = su xx + *u% 
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6.4.1 Deformation 

First of all we consider the term 



diXf \ d\X^ ] . By using the formula (2.20) we get 



-(2) 



du s (x) 

d Ux?A 



XV QS 



Ux?) 

V / : 



d (diX, 



(2) 
2 



\ / zv 



d 



du s (y) 

Ux?) 

\ / : 



+ 



yz 



du s (x) X \ 2 J: 



8u s (z) \ / zy UU-yXj \ / xy 



+ 



Let us focus our attention on the first term. By straightforward calculation we get 

d fax®) 

=l,...,3y'=l,...,6-i 



s 

<9w s (x) 



(6.17) 



(6.18) 
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where 




The other terms in (6.17) have the same form. The only difference is that the variables x, y and z 
play a different role. Therefore we can apply the usual tricks and write an expression containing 
only terms with 5 <yl \x — y)8^\x — z). For example we can write 



34 



d s y Ux?) = - E Mv)* (i) (* - v)s®(x -z) = 



du s (y) 

"J 

ij 

= -(-iy +3 J2 b M d ^ - yW x - z )) 

ij 



fc=0 



(-ir +i ^^.(y) E u ) * (fc) (* - v)* '-*- fc) (* - *) = 



(-!) m EE( I )E( k i |M^ N (»-y)^-^) 



The final result is 



<9s 



ij fc=0 v 7 Z=0 



16 ( S ^L\ s 1 A * 5 + 40 (s|^ A S 1 + 403, (s|^ i 1 A S 1 + 



+m ' I S SI J * 2 A 153 + 329= ( S g A (S 3 + 8C4? (,| ) f ■ A «» 
The calculation of the term d2X^ can be done as above. In this case we have 



doX. 



(4) 



2^-2 



xy 



-jd x X + cio + c 2 i + c 32 + C43 + c 54 ) 5(x - y) 



(X + c 00 + c\\ + C22 + c 33 + c 44 ) <5 (1) (:e - y) + (coi + C12 + c 23 + c 34 ) 5 {2 \x - y) 
+ (c 02 + C13 + c 24 ) 5 (3) (a; - y) + (c 03 + c 14 ) c5^(x - y) + c 04 5 {5 \x - y) 



where X { 2 A) = X = 



squ xxxx + s\u x u xxx + S2(u xx ) 2 + szu x u xx + S4U* and the coefficients c^- have the 
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same expression in terms of X^f 1 that in the previous case, except for the new coefficients 
c 54 = -ud* f-^U - l -^u - l -u 9 4 

(^54 au x \ f) I o a u xxxxx <-> u x u x I o 

\OU X xxx J A VU'xxxx A \V" J xxo 

c 4 4 = -5uo x ( - ) - 3- Uxxxx - 2u x d x 



r\ I £k ^XXXX ~^X^x \ f) 

L/ ^XXXX / XXXX \ 

c 34 = -IOmcL - 7- u X xx - 3u x o x 



r\ I • r\ ^XXX ^^X^x \ O 

yj f-t/xxXX / ^ ^XXXX \ ^ ^XXXX 

c 24 = -10ud x I I - 8- u xx - 2u x d x I 

V OU X xxx J OUxxxx V O^xxxx 

/ 31 \ <9X 
C14 = -5ud x I 7; I - 5- u x 

\ OU xxxx J OU xxxx 

c 04 = a 

oU xxxx 



Consequently we can write 



(d 2 xfA =<T,b k 6M(x-y) 

\ / xy ~ — 

with degipk) = 5 — k and 



xy 

k=0 



bo = a u X xxxx + aiUxUxxxx + a 2 u X xUxxx + a 3 u 2 x u xxx + a A u x (u xx ) 2 + a b u x u xx + a 6 u x = 

ds \ ( ds ds 2 dsx d 2 s \ 

2us 1 - 2us 2 - 2u— Uxxxxx + si - -a s 2 - 6m — — + 6m— 6m— — u x u xxxx 

du J \ du du du du 1 J 

f in 3s 2 9si <9 2 s \ ( dsi ds 2 d 2 s ds 3 
+ -10m— + 10m- 10m— — UxxUxxx + 2- 2— 2—— + 12ms 4 - 4m- — h 

\ <7M <7M 0M 2 / \ OU OU OU OU 

d 2 s 2 d 2 Sl d 3 s \ 2 / 3 9s 2 3 9 Sl 39 2 So ^ 9/1 a ds 3 

~ QU ^ + ~ ) U * U ™ + [-2 ~du~ + 2 ft7 " 2 ft? + 24MS4 ~ 8M ft7 + 

&S2 &S! d 3 s \ \2,( R ds s d 2 s 2 d 2 Sl <9 3 s 
—- - + 15m—— — 15m - - tir Mtt I + 6S4 — 2— — — + 6—— — 6 — — ■ 



~ 7u ^~Y + lbu ^~Y ~ 15m ^T "^(Wxx) + 6s 4 - 2- — — + 3— T - 3— —-\ 

ou z ou z ou 6 J \ ou OU A ou z ou 6 

ds^ n d 2 s 3 n d 3 s 2 t _ tn (9 3 si _ tn d 4 s \ 3 Z' 3 <9s 4 l<9 2 s 3 § 1 9 3 si 



+24M ft7 " 8u ^ - 2u ^ + 10u ^ ~ 10u ^ u * u - + (2^7 - 2^ + 2^ + 



l d 4 s g 2 s 4 a 3 s 3 >9 4 si 5 
'2ft? + 3M ftI 2 "" M ^ + M ^~ M ft?y' % 
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bi = c u xxxx + c x u x u xxx + c 2 {u xx ) 2 + c 3 u 2 x u xx + c 4 u x = 

-2s - 7m^ + Qus l - Qus 2 ) u xxxx + ( -2s 2 - s 1 - 2^ - 2ms 3 - 12m^ + lQu^p-+ 
ou J \ ou ou ou 

-20m^^^ u x u xxx + ^-3s 2 - 2ms 3 - 8«t^ + 12m ~^ ~ 15m ^^) ^ Uxx ^ + (~ 5s3 + 

<9s 2 , 9 9si <9 2 s <9s 3 R <9 2 s 2 . OA <9 2 si <9 3 s \ 2 

-2— + -— 6— — + 24ms 4 - 13m- 6m— — + 24m—— - 30m— — u x u xx + 

OU 2 <7M OM^ <7M ou 2 ou 2 OU A J 

ds 3 3d 2 s 1 <9 3 s . a ds A d 2 s 3 d 3 s 1 d 4 s . 4 

ou 2 air ou 6 ou ou 1 ou 6 ou^ 



x 



h = d u xxx + diu x u xx + d 2 u\ = ^-6ms 2 + 6us 1 - 10m ^ ~ 7s o ) u xxx + 

/ ds ds 2 ds 1 d 2 s \ 

+ -6s 2 - 3si - 3— 6ms 3 - 6m— + 18m— 30m— — u x u xx + 

V ou ou ou Ou 1 



3 ds 1 <9 2 s <9s 3 <9 2 si <9 3 s 1 

-3s 3 + -— 3^—^ — 3m— h 6m ^— - — 10m - „ ) m 



3 

,,,, v , , ( I V , I ( 1 " 

2 du du 2 du du 2 du 3 



h iH)a,,v f -8s - 4ms 2 + 4msi - 10m^ ] //,,,. 4 



2 



ds A ds 1 d s , 

-3si - 2— 2ms 3 + 4m— 10m—— u 

ou ou ou 2 J 

h = qu x = f-5s - 5u^-^j u x 

65 = k = — 2ms 

We are almost able to write the equation 

diX^diXf)] = (6.19) 



d\d 2 X^ ^ 



2 

in terms of b^. In fact to obtain the term did 2 X^ it is sufficient to add the following terms to 
the formula we have found before: 

\^5^{x - y)6W(x -z) + l^ 5 \x - y)6W(x - z) + 

Zi (J LI Zt (J LI 

+\^-S^\x - y)6W(z -z) + \^^\x - y)S^(x - z) 

— C/ ^ XXX xxxx 

+lJ*-8(x-y)6V>(x-z) + ... 
z ou xxxxx 

and to rewrite the terms in the usual way. 
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Defining the functions /, g, h in the following way: 
db 5 _ <% 4 + db 2 <%i 



(^H/n^ ^^ r ^ j xxx ^-^^-^ xxxx 

du xxx du xxx du x du \du xxxx I \du xxx J \du x I \du 



du x x \du xx ) x \du xxx ) x \du xxxx J du x \du) x \du 



x 1 du 



we have 



d x d 2 X = 2/5 1 A5 5 + 5f5 2 A5 4 + 2g5 A 5 5 + (5f x + 5^)5 X A 5 4 + (6/* + 2#)5 2 A 5 3 + 5^<J A 5 4 
+ {4f xx + %g x )5 l A 5 s + (2h + Ag xx )5 A 5 3 + (f xxx + 3h + Zg xx )5 l A 5 2 + {3h x + g xxx )8 A 5 2 + 
+h xx 5 A 5 1 



Then the equation d\d 2 X^ — | d\X^ ) ', g^X^ = is eqivalent to the conditions 

/-4^ = 
aw 

2 = 
/i = 

The first equation connects the second order deformation with the fourth order deformation: 

d s 

s = -2s— (6.20) 

The second and the third equation give the conditions 

ds 



(2) A V( 2 ) 



Si - S 2 = -TT - (6.21) 
du 

ds 3 d 2 s 2 ( . 

354 - ^ + = ° (6 ' 22) 



6.4.2 Trivial deformations 

In this case trivial deformations can be written as d\A + d 2 B with 



A= / (A (u)u xxx + Ai(u)u x u xx + A 2 (u)u 3 x )dx (6.23) 
B= (B (u)u xxx + Bi(u)u x u xx + B 2 (u)u 3 x )dx (6.24) 
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A straightforward calculation shows that, for trivial deformations, the coefficients Sj are 

s = (6.25) 

Sl = S2 = 3A + 3uB (6.26) 

OA „ dB 3~ ,„ rt „, 

S3 = 6 ^ + 6 "a7 + 2 B (6 ' 27 > 

with 

dA 1 d 2 A 



A = 2A 2 7T- 4 

5 = 2B 2 - ^ 4 



du du 2 
dB l d 2 B 



du du 2 

Now we are able to prove that P 2 4 ^ is trivial if and only if s = 0. 
Proof 

Trivial deformations satisfy equations 

s = 



Sl = s 2 
c 

du ' du 2 



ds 3 d 2 s 2 

3s 4 h -t^t = 



In fact 



,2 



du 2 V / 2 du du \ du V / 2 / <9w : 

Conversely, we can always write a deformation X (with s = 0) as X = d\A + d 2 B. It is sufficient 
to choose A and .£> in such way that the following equations hold: 

5 Ads 2 2 

B = -— Ss 

9 du 9 3 

A = -s 2 -B 
3 

6.4.3 Deformations: explicit form 

By using the formula (2.19) 
d x X = 

( d 
d 

~ (9 s J'(x - y)) du(s) ,J s u xxxx + s x u x u xxx + s 2 «L + s z u t) 



8 to avoid any confusion with the other functions s defined in this paper, in the theorem 4 we have substituted 
sq for s. 
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and conditions (6.20), (6.21) and (6.22) it is easy to get the formula 



d,X^ = -2s 5^(x -y)-5 (d x s ) 5^(x - y) - 10 (d 2 x s ) 5^(x - y) + 

-10 (figso) S (2 \x - y) - 3 (%s ) 5 (1) {x - y) + 2^5^ (x -y) + 3 (d x w) 5 {2 \x - y) + 

+ {d 2 x w) 6W(x-y) 

whith 

w = w u xx + w^l = 2^-u xx + Wiii 2 , (6.29) 

where w\ is an arbitrary function of u and so is related to the function s appearing in second order 
deformation by the equation 

ds 

s = -2s— (6.30) 



7 Quasi-triviality 

Definition 16 The group of transformations 

k Fk\ u ) U X) u xx) ■■■) fn i\ 

u -> u = 2_^e — (7.1) 

Ljr k [U, U X) U XX) ...) 

where F k , G k £ A, deg(F k ) — deg(G k ) = k and ^ ^ is called quasi-Miura group. 

A deformation is called quasi-trivial if it can be eliminated by the action of the quasi-Miura group. 
We have seen that, in general, the second and fourth order deformations are not trivial. In this 
section we show that they are quasi-trivial. 



7.1 Second order 

Theorem 6 If = d\A + d 2 B, where A and B are local functional whose densities are ratios 
of differential polynomials, then the secon order deformation P\ = P^ — X(P 2 °^ + ^d\X^) is 
trivial. 



Proof 

Xf ] = d x A + d 2 B =>■ P 2 (2) = -d 2 d l B 
This implies 

P 2 (2) = Lie x P^ 
Lie x Pi ] = 

with 

X = -d x B (7.2) 
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Theorem 7 All second order deformations are quasi-trivial. 



Proof 

If we choose 



A 
B 



I (a(u)—^-)dx 
(b(u)—^-)dx 



By using formula (2.18) we get 



d x A + d 2 B = 



' du 2 



+ 2u 



d 2 b 1 db 
du^ + 2du~ 1 



da db\ u 2 
tt- + u— 1 ' 

du 



du ) u 2 x 



+ 



(da db\ u xa 
\du du J u a 



+ 



d 3 a d 3 b 1 d 2 b 
+ u— + 



du 3 



du 3 2 du 2 



u , 



If we put 



Then 



a(u) = —ub(u) + 

J u 



b{u)du 



db 

du du 



da 

ht 
d^a 

du 2 
d^a 

du 3 



db 



d 2 b 



u- 



= -2 



du du 2 
d 2 b d 3 b 



du 2 



u 



and these equations imply 



d x A + d 2 B = -- I -t^Uxx + 



du 3 

d 2 b 
du 2 



u„ 



that is equal to X 2 



(2) 



s(u)u xx + ^u 2 x if we choose 



b(u) 



s(u)du 



UQ 



7.2 Fourth order 



(7.3) 
(7.4) 



(7.5) 



(7.6) 



(7.7) 



Lemma 1 A fourth order deformation is quasi-trivial if and only if there exists a vector field Y 
such that 



Lie Y Pl 0) = (7. 

where X is the vector field (7.2). 



P 2 (4) - \hie\P^ = LieyPW (7.9) 
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Proof 

The reduction of the fourth order deformation to the form P^ — XP 2 ^ can be achieved in two 
steps. 

In the first step one kills the second order part of the deformation P 2 ^ with the quasi-Miura 
transformation generated by the vector field (7.2): 

r* - > (* + ^ + + o ( „) . p« - » (* + , (tf> - i^f + 0( ,,) ) ^ 

In the second step one kills the fourth order part of the deformation P^ — ]^Lie\P 2 ^ with the 
quasi-Miura transformation generated by the vector field Y . 

Theorem 8 All fourth order deformations are quasi-trivial. 

To prove the theorem we need the following 
Lemma 2 If = d x A + d 2 B - l[d±B, d 2 B] where 

B= [ b(u)—dx 

and b(u) is given by the expression (7.7), then the fourth order deformation 

is trivial. 
Proof 

P 2 {4) = rfiXf } = d x d 2 B - d x \\d x B, d 2 B] 
By using the graded Jacobi identity we get 

P 2 (4) - \ue\P 2 ^ = Lie Y pW (7.10) 

where X = —d\B and Y = —d\B. Moreover 

Lie Y P ( i ] = 

Then to prove the theorem it is sufficient to show that the vector field X^f 1 = s u xxxx + S\U x u xxx + 
S 2 U XX + s 3 u 1 u xx + s 4 u t 9 can De written in the form 

X^ = d x A + d 2 B - \[d x B, d 2 B] 

9 We recall that the coefficients of the vector field satisfy the condition (6.20), (6.21) and (6.22) 
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We start calculating the term —^[diB, d 2 B]. In this case the Schouten bracket coincides with the 
commutator of the vector fields 



dB = 2—u — ^ I dbUxxx | d 3 ft u 2 

du 2 xx du u 2 du u x du 3 x 

and 

1 fdb d 2 b 



diB = ud lB + -\^u xx +~—ul 



The result is 




27 / d 2 b \ 2 45 db d 3 b 



4 <9w du 3 



+ 



, 45 d 2 bd 3 b 21 db d 4 ^ 

, . 9_&bd*b Sffibdb 3 /^6\ 2 \ 9^ /96\ 2 
* l 4<9w 2 ,9w 3 2 9m 5 9m 4V<9mV y 2 u\ \du) 

g}4x ( 96 V + 2 W / 9&y _ ^m^x / 27 96 d^A 
u s J u x \du ) u x \2 dudu 2 ) 

u\ x ( 33 db^b\ ,//;,,.//,,,,,. / d6\ 2 _ Uxx^xxxx ( db\ 2 



4 aw aw 2 / u% \du J uj. \du J 



The problem now is to guess the form of the local functionals A and B. We know that the degree 
of A and i? is equal to 3 but all ratios of differential polynomials of the form 

P(u, u x , u xx , ...) 



Q\U"> "i) u xx , ...) 

where deg(P) — deg(Q) = 3, are allowed. 

The form of the coefficients in (7.11) suggests that A and B contain only terms of the form 

P(u,u x ,u xx ...) 

K 

for some integer i. 
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Let us try with the local functionals 



f fii / \^xx i i in/ \^xx , ; nl x^xx^xxx , A , \U XXXX , lr , \U XX U XXX , \U xa 

a= [hliu)-^ + +h2(u)-^ + hZ(u) h M(u) h h5(u)^— h/iom)-^ 

J s i V u x ul u 2 x u x u% '/;;. 

. r \U XX U XXXX , ,„/ lUnira , -, A / \U XXX U XXXX . <.U \ , 

+h7(u) h h8{u)^± L + h%u) - h hlOiu) -. h hlliu) , eta; 

«s «x M x < «| / 



• 1 ( // ) — + +£;2(u)% + k3(u)^^ + kA{u)^^ + fc5(u)^^ + k6(u)^f+ 
u x u x u x u x u x u x 



+k7{u) UxxU l xxx + k8{u)^f + A;9( M )^= + ' , ' + fcll( u )%^ > ) da: 

M X U I U x U x U x J 

In order to determine the exact form of the coefficients hi and ki one has just to compare did + d 2 b 
with \[d l B,d 2 B\. 

The first step is to kill the coefficients of rational terms in d\d + d 2 b that don't appear in 
— \[d\B, d 2 B] and the second step is to eliminate the remaining rational terms. 
The calculations are very long, we give here only the final result: 10 

hi(u) = —uki(u), i — 1, 11 

fcl(«) = -^*8(u) - ^*2(u) 

k3(u) = Akl(u) 
M(u) = 
fc5(u) = -6k8(u) 
k6(u) = Ak8(u) 
k7(u) = k8{u) 

fc9(u) = fclO(u) = fcll(u) = 
, . , dk8 7 (db\ 2 

Whith this choice and taking into account (7.7), by using formula (2.18) we get 



--[diB, d 2 B] + did + d 2 b 



-2s ^ U xxxx +| rt -5* n ( n ) I ^x^a 




ds f 7 d 2 s 13 /9s 

(9-u 



,'5 <9 3 s ds<9 2 s\ 2 /l/<9 2 sY dsd 3 s 2 <9 4 s 



/ _|_ / , / 2 , / 2 _|_ / 4 

SqU xxxx -+- S^U X U XXX ■+- S2U XX -+- S^U X U XX ■+- <S4li a 





10 In fact there is a certain freedom in the choice of coefficients /i, and fcjj for example it is not necessary to put 
A;4(u) equal to 0. 
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First of all we observe that s' is equal to so- Moreover the coefficients s^ satisfy equations (6.21) 
and (6.22): 



d 



Sl 82 du \ & du 



ds 



4 du du 2 



d 2 s\ 2 ds d 2 s d 4 s 5 d 4 s 

3—— + 2s— - -s— + 



du 2 



du du 2 



7dsd 3 s ( d 2 s\ 2 16 ds d 3 s fd 2 s^~ 



3 du du 3 



du 2 



3 du du 3 



du 4 3 du 4 
1 d 4 s 



du 2 



3 du 4 



Then the vector field 
1 



X - (--[^B, d 2 B\ + d x d + d 2 b) 



1 \.A 

X 



= (so - s' )u xxxx + (si - s[)u x u xxx + (s 2 - s' 2 )u 2 xx + (s 3 - s' 3 )u 2 x u xx + (s 4 - s' A )u 
is trivial. Indeed: 

so - s' = 



Moreover the equations 



si - s 2 
si - s 2 



ds 

du 
ds 

du 



imply 



Finally 



Si — Si — s 2 — s 2 



d d 2 
3(s 4 - s' 4 ) - — (s 3 - s' 3 ) + — (s 2 - s' 2 ) 



3s 4 - 77- s 3 



d 2 



du du 2 



s 2 - [3s'a- 



d_ , cP_, 
du S3+ du 2 " 2 



This means that 



that is 



X - (--[diS, d 2 5] + tZia + rf 2 &) = d ± a + d 2 b 



X = --[diB, d 2 B] + d ± (a + a) + d 2 (b + 6) 
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8 Conclusions 



In this paper we have studied the problem of classification of deformations of bihamiltonian struc- 
ture of hydrodynamic type. 

The main result of the paper is that, up to the fourth order, any deformation can be reduced to 
the form (3.15) depending only on a functional parameter s(u). 

These deformations give rise to an infinite hierarchy of "almost-commuting" hamiltonian equa- 
tions. 

In this paper we started studying numerically one of the equations of the deformed hierarchy 
corresponding to second order deformations. 

The results obtained indicate the existence of the analogue of 2-soliton solutions at least for small 
times and for small amplitudes, but a deeper analysis is still necessary. 

We have seen that one of the consequence of classification theorem is that deformations of the 
Magri bihamiltonian structure are trivial up to fourth order. This fact suggests that in our classi- 
fication problem the Magri bihamiltonian structure is a stable object according to the definition 
of Arnold. 11 It would be interesting to investigate if this situation is more general, that is, 
if completely integrable systems of bihamiltonian type are the stable objects of a corresponding 
classification problem. For example it is possible to apply the same techniques used in this paper 
to the "symplectic" case where the leading term in the parameter e has the form 

P li2 = h y(x - y) (8.1) 

with 



h\ j = 



/o 


1 





.. 


. 


o\ 


-1 








.. 


. 














1 .. 


. 











-1 


.. 


. 














.. 


. 


1 


\o 








.. 


. -1 


0/ 



and 

11 A classification problem can be thought as a decomposition of a space of objects in equivalence classes (see 
[2]); in our case the objects are deformations; two deformations are equivalent if and only if they can be reduced 
to the same form by the action of Miura group. An object is called stable if a sufficiently small "neighbourhood" 
of this object contains only objects of the same class. 
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Also in this case the the cohomology groups H 1 and H 2 associated to the differentials dp 1 and dp 2 
are trivial (see [7]) and then the non trivial deformations are classified by the group 



Finally we observe that the result about quasi-triviality suggests that this property is a conse- 
quence of the definition of deformation and it is not an additional constraint. 
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